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ABSTRACT 



' Context. In coded mask techniques, reconstructed sky images are pseudo-images: they are maps of the correlation between the image recorded 
, on a detector and an array derived from the coded mask pattern. 
Aims. The INTEGRAL/IBIS telescope provides images where the flux of each detected source is given by the height of the local peak in the 
correlation map. As such, it cannot provide an estimate of the flux of an extended source. What is needed is intensity sky images giving the 
flux per solide angle as typically done at other wavelengths. 
O l' Methods. In this paper, we present the response of the INTEGRAL IBIS/ISGRI coded mask instrument to extended sources. We develop a 
general method based on analytical calculations in order to measure the intensity and the associated error of any celestial source and validated 
with Monte-Carlo simulations. 

Results. We find that the sensitivity degrades almost linearly with the source extent. Analytical formulae are given as well as an easy-to-use 



C/2 

, recipe for the INTEGRAL user. We check this method on IBIS/ISGRI data but these results are general and applicable to any coded mask 
• • 1 telescope. 
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1. Introduction 

The usual techniques of imaging the sky by focusing the radia- 
tion become ineffective above ~ 15 keV due to technological 
constraints. Future observatories (e.g. Constellation-X/HXT, 
Hailey et al. 2004 ; Simbol-X, Ferrando et al. 2005) will push 
this limit up to 60-70 keV wi th the use of multi-layer graz- 
ing optics jRamsev et al. 20021) . Collimators using the standard 
on/off techniques are limited by the source density at low en- 
ergies (< 50 keV) and by the internal background variability 
at higher energies. Both problems are addressed with coded- 
aperture techniques, which have been employed successfully 
in high-energy astronomy in several previous mission s such 
as the GRANAT/SIGMA instrume nt ipaulet al. IQQll) or the 
BeppoSAX Wide Field Cameras lla geret al. 1997h . In such 
telescopes, source radiation is spatially modulated by a mask 
consisting of an array of opaque and transparent elements and 
recorded by a position sensitive detector, allowing simultane- 
ous measurement of source plus background in the detector 
area corresponding to the mask holes and background only 
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in that corresponding to the opaque elements. For each point 
source in the field of view the two-dimensional distribution of 
events projected onto the detector reproduces a unique shadow 
of the mask pattern and the recorded shadowgram is the sum 
of many such distributions. The popular way to reconstruct 
the sky image is based on a correlation procedure between 
the recorded image and a decoding array derived from the 
spatial characteristics of the mask pattern. One can describe 
the usual properties of telescopes for these systems as fol- 
lows: the angular resolution, defined by the angle subtended 
by the smallest mask element seen from the detector, and the 
field of view (FOV) divided in two parts: the Fully Coded 
Field of View (FCFOV) representing the sky region where the 
recorded source radiation is fully modulated by the mask and 
the Partially Coded Field of View (FCFOV) the region where a 
source projects only a part of the mask pattern on the detector. 

Many mask patterns have been designed and employed 
in the past for optimizing the imaging quality of point-like 
sources in the high-energy domain. Recently, Schafer (2004) 
has presented a novel method for imaging extended sources 
by constructing mask patterns with a pre-defined Point Spread 
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Function (PSF). In this paper, we have developed a method 
to reconstruct the intensity and estimate the corresponding er- 
ror of extended sources seen through a given coded mask tele- 
scope. This work is motivated by the fact that several sources 
of interest in y-ray astronomy must be considered as extended 
in order to properly derive their intensity such as supernova 
remnants (hereafter, SNRs): the Tycho SNR (with an appar- 
ent diameter of ~ 8'), SN 1006 (0 ~ 30'), G 347.3-0.5 
(0 ~ 1°) or RX J0852-4622 {0 ~ 2°), clusters of galaxies {e.g. 
Coma Cluster, ~ 1°), and diffuse interstellar emission from 
various high-energy processes. We have applied this method 
to the IBIS/ISGRI instrument onboard the INTEGRAL satellite 
but the results presented here could be easily adapted to any 
coded mask telescope. We briefly describe the basic properties 
of IBIS and the principles of the standard data analysis we used 
for simulating celestial extended sources in section|2j consider- 
ing the simplest case of a uniform disk. We present our results 
concerning the IBIS response for this source geometry over a 
large range of radii (section ITlT l as well as in an astrophysi- 
cal case: the supernova remnant SN 1006 ( section l3T2t . These 
results allow us to find a general method to reconstruct the 
flux and the associated error of an extended source presented 
in sections WA\ and |431 respectively. We also describe the tests 
we performed on the Crab Nebula and compare the results of 
this method to those of the standard analysis (section l4!2t . In 
a joint paper, we have successfully ap plied this method on the 
IBIS/ISGRI data of the Coma Cluster iRenaud et al. 20061) . the 
first extended source detected with IBIS. Finally, we give prac- 
tical tips on the use of this method for any interested observer. 

2. Imaging sources with a coded mask telescope 

2.1. General imaging procedure 

In coded aperture imaging systems, a mask consists in an array 
M of 1 (transparent) and (opaque) elements. The detector ar- 
ray D (the real image) is simply the convolution of the sky im- 
age S with M plus an unmodulated background array B. If there 
exists an inverse correlation array G such that (M*G),j = Sij, 
then the reconstructed sky image is given by the following op- 
eration jFenimore & Cannon 198ll) : 

S' ^D-kG^S i^M-kG + BirG^S +8*0 (1) 

In the case where the mask M is derived from a cyclic repli- 
cation of the same basic pattern and the background B is given 
by a flat array then the term B*G is constant and can be re- 
moved. The quality of the objec t reconstruction th erefore de- 
pends on the choice of M and G jCaroU etal. 19871) . Fenimore 
& Cannon (1978) have found mask patterns with such prop- 
erties, including the Modified Uniformly Redundant Arrays 
(hereafter, MURAs). These coded masks are nearly-optimum 
and possess a correlation inverse matrix by setting G - 2M - 1 . 
For those systems the resulting sky image of a single point-like 
source in the FCFOV will have a main peak at the source po- 
sition of roughly the size of one projected mask element, with 
flat sidelobes in the FCFOV (for a perfect detector) and 8 main 
source ghosts in the PCFOV. 



2.2. IBIS properties 

We consider the IBIS telescope dUbertini et al. 2003h 
launched onboard th e ESA y-ray space mission INTEGRAL 
jWinkler et al. 2003^ on October 2002. The IBIS coded mask 
imagi ng system includes a replicated MURA (see section IZTl 
mask jGottesman & Fenimore 198^ of tungsten element s and 
2 pixellated y-ray detectors: ISGRI jLebrun et al. 20031) . the 
low-energy camera (15 keV - 1 MeV) and PICsIT , operating 
between 175 keV and 10 MeV jLabanti et al. 20031) . We focus 
here on IBIS/ISGRI but the results can be adapted to any coded 
mask imaging system. The physical characteristics of the IBIS 
telescope define a FCFOV of 8.3°x 8.3° and a total FOV of 
29° X 29° at zero sensitivity. The nominal angular resolution 
is 12' (FWHM) and standard ISGRI images are sampled every 
~ 5' , which is the angular size of detector pixels seen from 
the mask, allowing fine imaging and positioning of detected 
celestial sources. The decoding array G is oversampled at 
the pixel scale in order to optimize the signal to noise ratio 
of point-like sources in the reconstructed image. One may 
refer to the papers of Gros et al. (2003) and Goldwurm 
et al. (2003) for details on the shape of the System Point 
Spread Function (hereafter, SPSF) and the full description 
of the implemented algorithm for IBIS data, respectively. 
Goldwurm et al. ( 2003) show a variation of finely balanced 
cross correlation jFenimore & Cannon 198lh generalized to 
the total (FC+PC)FOV. 

2.3. Imaging extended sources 

In the following, any extended source will be considered as a 
large number of unresolved points. The principle of our method 
for imaging extended sources is as follows: first, we have cal- 
culated the fraction of pixel area covered by a given projected 
mask element (called Pixel Illumation Factor, PIF) correspond- 
ing to each point-like source given its relative sky position. 
These calculations take into account all the instrumental fea- 
tures of IBIS such as efficiency, dead zones between pixels and 
detector modules of ISGRI, mask thickness and transparen- 
cies of all the intervening materials in the ISGRI energy range 
(15keV - 1 MeV). Second, we have performed a weighted sum 
of each of the point-like source contributions to obtain the shad- 
owgram of any extended sky source according to the following 
equation: 

N 

Du^Y,f«''P^^u (2) 

1V=I 

where is the image recorded on ISGRI, N refers to the 
number of unresolved points forming the extended source, and 
f„. represents the flux of each point source (w) with 2 fw - F, 
the total flux of the extended source. 

In order to illustrate the blurring effect due to the source 
extent. Figure [0 shows the ISGRI shadowgram of PIF values 
corresponding to an on-axis point-like source and that of a 
uniform disk of ~ 20' radius. The overall pattern of the PIF 
shadowgram has a very little dependence on energy and only 
the relative amplitude of fluctuations depends on the energy. 
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Fig. 1. The recorded ISGRI shadowgram of the PIF values for 
an on-axis point-Hke source (left) and that of a uniform disk of 
~ 20' radius (right) without background. Dead zones between 
modules and mask elements are clearly visible. 



Finally, in order to reconstruct the sky images, we have used 
the standa rd deconvolution inclii ded in the Off-Line Scientific 
Analysis dOoldwurm et al. 2Q03I) . version 5.0. We decided to 
perform these simulations in connection with what an observer 
can expect by analyzing IBIS/ISGRI data of any celestial re- 
gion. We have chosen to simulate a disk as discussed in sec- 
tion^ assuming that each individual point source contribution 
forming the extended source emit the same flux (i.e. the factor 
fw in the equation |2l is constant, each extended source is uni- 
form). We have distinguished two different relative positions in 
the sky: an on-axis source and one in the PCFOV, 9° away from 
the telescope axis, i.e. where the IBIS sensitivity is reduced by 
a factor of almost 2. 



3. Sensitivity to extended sources 

Simulations were performed for several radii, from {i.e. , a 
point-like source) to ~ 1° keeping constant the total flux F 
equal to YI^^.^i f» - should notice that in the standard re- 
constructed images, the flux of point-like sources in the FOV 
is given by the peak of their associated SPSF after fitting with 
a two-dimensional gaussian (see Gros et al. 2003). Therefore, 
for each reconstructed sky image, we have measured the max- 
imal pixel flux within the apparent diameter of each extended 
source (hereafter, f„,ax) to estimate the IBIS/ISGRI response to 
the global peak reconstruction. 

3.1. Uniform disk 

This case could relate to the Crab-like SNRs (plerions or 
pulsar-wind nebulae) containing pulsars that inject a relativis- 
tic wind of synchrotron-emitting plasma into the surrounding 
medium. These pulsar-wind nebulae can be seen in first approx- 
imation as uniform disks. Figure |2lpresents fmax as a function 
of the radius for the two different source positions in the IBIS 
FOV. 

The two curves present different behaviors; below a radius 
of ~ 8', the extent of the source is still comparable to the IBIS 
angular resolution and the flux loss is less than 30 %. For more 
extended sources, the SPSF differs from that of a point-like 
source and the flux decreases as R"^ (and R ' in the case of 
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Fig. 2. The relative reconstructed peak flux as a function of 
the radius of the uniform disk for two different positions: on- 
axis (triangles) and 9° away from the pointing (labeled off-axis, 
squares). The two curves have the same behavior and present 
two main features: below ~ 8' radius, the loss in flux is limited 
whereas above this value, the flux decreases as R"^. The dotted 
lines give the expected behavior based on the equation|3 



a uniform ring) due to the simple flux dilution. One can easily 
show that this is well described by the equation: 



fm 



x(l 



-K-/20-2, 



'0 



(3) 



where Cspsf is the width of the SPSF and R the source ra- 
dius. The difference between the on-axis and off-axis curves is 
coming from the fact that the width of the IBIS SPSF is not 
constant over the relative positions in the FOV: we have fitted 
the SPSF for the two positions with a bi-dimensional gaussian 
and found that the width in the on-axis case is close to the nom- 
inal value, ~ 14' (FWHM) and is around 17.5' in the off-axis 
one. These results are similar to those of Gros et al. (2003, Fig. 
2) and explain what we observe in Figure |2 the reconstructed 
peak flux of an extended source depends less sensitively to the 
source size when the width of the SPSF is larger We will dis- 
cuss in details the characteristics of the IBIS SPSF in section 

1331 



3.2. An astropliysical case: SN 1006 

SN 1006 is an historical Galactic SNR of ~ 30' diameter, ob- 
served with INTEGRAL for ~ 1 Msec in the first two years of 
the mission. Reynolds (1999) has modelled the predicted hard 
X-ray emission of this SNR with two components: the syn- 
chrotron radiation, concentrated in two bright limbs and based 
on the calculations of Reynolds (1998), and a non-thermal 
bremsstrahlung component more symmetrically distributed. At 
energies below about lOOkeV, synchrotron emission will dom- 
inate. Figure |3]presents four images of this SNR, with the ex- 
pected morphology of SN 1006 at 28 keV, the sampled image 
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in the ISGRJ sky pixels, and the resuhs of our simulations. In 
this example, SN 1006 is located in the FCFOV. 




Fig. 3. (a) The expected image at 28 keV (in units of 10" 
ergs cm"^ s"' Hz"' sr"') based on simulations performed by 
Reynolds (1999). (b) The sampled sky image in the ISGRI sky 
pixels size (in units of ct s"') . (c) The corresponding shadow- 
gram of the PIF values and (d) the reconstructed image after 
the standard deconvolution (in units of ct s"'). 

The results of the analysis of the JEM-X and IBIS/ISGRI 
data of SN 1006 are presented in Kalemci et al. (2006). The 
source was not detected with ISGRI. Using the same method 
as in section im it was possible to extract an upper limit on 
the hard X-ray emission of SN 1006. Since the expected mor- 
phology is almost symmetric and the emission is concentrated 
in two extended limbs, the ratio between the maximal value in 
the reconstructed image at the position of each limb and their 
respective global flux provides a correction factor. In Kalemci 
et al. (2006), the "standard" upper limits (obtained as if the 
SN 1006 is an unresolved source) were divided by this correc- 
tion factor of ~ 0.7 to take the flux dilution into account. Note 
that a relative reconstructed peak flux of 0.7 is close to what we 
obtained for an extended source of ~ 10' radius (see Fig.|2j- 



biased estimation of the flux of an extended source. For proper 
flux estimation, one needs to convert the reconstructed standard 
images into images of intensity (i.e. flux per solid angle or sky 
pixel), as is typical in astronomy. In that case, the global flux of 
any source (point-like or extended) would be given by integrat- 
ing intensities over the source extent. The method is to divide 
standard flux images by the SPSF effective area (the solide an- 
gle subtended by the SPSF on the celestial sphere, i.e. the ratio 
between its integral and the peak value) but the main difficulty 
is that the shape of the IBIS SPSF varies within the FOV, so it 
is necessary to construct a map of its effective area (hereafter, 
A). 

We have calculated point-like sources each ~ 15' (3 ISGRJ 
sky pixels sampling) in the first quarter of the IBIS FOV with 
a fixed flux and fitted each corresponding reconstructed image 
at the source position by a 2D gaussian. We assume that these 
images must have a central symmetry due to that of the MURA 
coded mask and then we have projected each quarter of im- 
ages where we performed the simulations on the three others 
according to the central symmetry. Figure |3 presents the map 
of A of the IBIS SPSF and the histogram of its values. The map 
is almost flat in the FCFOV and the values globally increase in 
the PCFOV. Its histogram shows an appreciable dispersion (cr 
~ 1.3 sky pixel") with a mean value of 8.9 sky pixel". To obtain 
reconstructed sky images in intensity, one can then apply the 
following transformation: 



I's-Yj^L, (5) 
n 

where S'osa is the reconstructed flux image obtained with 
the standard deconvolution. The legitimacy of such a transfor- 
mation is shown in Appendix lA. 1 1 and is based on the sole as- 
sumption that A does not change on scales smaller than the 
SPSF width. Thus, for extended sources simulated as described 
in section|2l we have built these intensity images /,'„,. For each 
of them we have integrated over the source extent as seen by 
IBIS (Q, i.e. the physical size of the source convolved with the 
IBIS SPSF ~ source radius plus 3 xcTspsf) and we have checked 
that this final value corresponds to the input global flux of our 
simulations with a dispersion of less than 10 %. 



4. General method 

We have shown that any coded mask imager such as IBIS works 
with extended sources just as predicted based on its SPSF prop- 
erties: the reconstructed peak flux follows the simple flux di- 
lution over the large number of unresolved points constituting 
an extended source. Moreover, as previously presented by Gros 
et al. (2003), the width of the IBIS SPSF varies within the FOV. 

4.1. Reconstructed flux of extended sources 

As described in section 13 the flux of a point-like source in a 
reconstructed IBIS image is given by the peak of the SPSF. 
Therefore, as shown in the previous sections, this will give a 



4.2. Tests on the Crab Nebula observations 

We have also tested this method on Crab Nebula observations. 
This source is considered as point-like for IBIS. We used data 
obtained during the revolutions 102, 170 and 300 containing 
164 science windows (individual pointings, hereafter sews) 
and applied the same procedure as described above. Figure 
|5l presents the relative positions of the Crab for the analyzed 
sews and the histograms of the fluxes found with our method 
and the standard one in the 20-30 keV energy range. Even if 
our method shows a larger dispersion in the reconstructed Crab 
fluxes {cr,„pii,„j ~ 3.1 s"' and ctosas.o ~ 2.3 s"'), the mean val- 
ues are compatible (<F,„e,;,o£/> ~ 78.9 s"' and <Fosa5.o> ~ 78.7 
s"'). Since the Crab Nebula is an unresolved source for IBIS, 
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Fig. 5. Left: Relative positions of the Crab Nebula for the 164 
analyzed pointings during the rev. 102, 170. Right: Histograms 
of the reconstructed fluxes of the Crab with the standard anal- 
ysis OSA 5.0 (in blue) and with our method (in black) 



the reconstructed flux with our method depends on the source 
location since A varies with the position in the field of view. 
However, the sensitivity of this dependence decreases with the 
source extent. 



4.3. Estimation of tlie associated errors 

Although the flux of any extended source can be measured with 
the method described above, the estimation of the correspond- 
ing errors is not simple. In fact, the reconstructed counts in the 
sky pixels within the SPSF of the IBIS telescope are correlated. 
Since we propose to sum the fluxes previously divided by A 
over the size of the source, the associated variance should be 
2 o-,'^„/A^,„ (o-'i^^ is the OSA flux variance) plus a covariance 
term. The latter depends on the position in the FOV and the 
source extent. 



We have developed the standard equations of the convolu- 
tion in Appendix IA.2l to find a simple analytical expression of 
the variance associated to the equations |3 and [S] such that: 



cr,; 



1 



F A 



a-,; 



spsf 



(6) 



(7) 



where cr^ . is the detector variance of the background count 
rate (in units of count s"-), E„„ is the imaging efficiency and 
< crj^^^ > is the OSA flux variance averaged over the region of 
pixel fluxes summation O. Nspsf refers to the number of SPSFs 
within f2. Equation0can be interpreted as follows: the width 
of the SPSF {(Tspsf) is the correlation length of the instrument 
beyond which sky pixels are uncorrected. Therefore, at these 
scales, the final variance associated with the reconstructed flux 
of an extended source by summing the pixel fluxes is simply 
given by the sum of the variances. This can be then approxi- 
mated to the mean variance times the number of SPSFs within 
Q. (see AppendixiBl. 

In order to check this analytical estimation of the associated 
error, we have performed Monte-Carlo simulations of extended 
sources seen through the IBIS coded mask as a sum of unre- 
solved points. Shadowgrams are given by equation|2]for which 
we have added a constant background term. The total flux of 
each extended source was fixed at 50 s ' and that of the back- 
ground at 10"' s"' . 500 shadowgrams are simulated for different 
source sizes (assumed to be uniform disks), from 0° (point- 
like) to 1° radius, following the equation: VD^/ x A^(0, 1) + Dki 
where D defines the ISGRI expected count rate shadowgram at 
the pixel (k,l). N(0,1) is the normal distribution. Each of these 
simulated shadowgrams was deconvolved using the standard 
method of the IBIS data analysis in order to obtain 500 re- 
constructed images in flux, variance and significance. We per- 
formed these calculations for three different source positions 
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in the FOV: one on-axis (Y=200.5, Z=200.5) and two in the 
PCFOV at ~ 9.5° (Y=300.5, Z=250.5, labelled "Off Axisl") 
and ~ 13.5° near the edge of the FOV (Y=120.5, Z=50.5, la- 
belled "Off Axis2"). These calculations allow us to control 
the statistics for any source extent and position in the FOV. 
We have reconstructed the source fluxes using the method de- 
scribed in section 14.11 and measured the width of each flux 
histogram, which corresponds to the error associated to this 
method. 




0.0 0.2 0.4 0.6 0.8 1.0 
Source Radius (deg) 

Fig. 6. The ratio between the measured error of our method and 
that predicted for a point-like source (ctq) as a function of the 
source radius for the three relative source positions. The solid 
line is given by the equation |6| and its large-scale behavior up 
to the total IBIS FOV is shown in the encapsulated plot. 

Figure|6lpresents the ratio between the measured error and 
that derived from the equation |6l for a point-like source (la- 
belled (To) as a function of the source radius for the three dis- 
tinct source positions. This ratio characterizes the loss in sen- 
sitivity due to the source extent. The error bars were calculated 
from the error of the fit performed on each reconstructed fluxes 
histogram to estimate the error of the method. For any source 
size and position, the dispersion is estimated to be less than 10 
%. 

5. Discussion 

We have shown that it is always possible to obtain an estimate 
of the flux of any extended source with coded mask imaging. 
This is an intrinsic property of these imaging systems where 
the source flux is obtained through a linear combination of 
detector-pixel count rates. We have derived analytical expres- 
sions for estimating the global flux and its associated error of 
any extended sky source which were tested by Monte-Carlo 
simulations. This work was done relating to what one could 
expect to observe in the y-ray domain through a MURA coded 



mask instrument such as IBIS onboard the INTEGRAL satel- 
lite. For IBIS, extended sources with a radius greater than ~ 
4.5° will be confused with their associated ghosts. One will 
need a specific method to clean the images but this is beyond 
the scope of this paper The advantage of this method, as long 
as the integration region contains all the source, is that the ob- 
server does not have to assume any source geometry during 
the analysis as is the case for methods using model fitting. Of 
course, in principle, the optimum signal to noise ratio is at- 
tained with the smallest integration region still containing all 
the source. However, due to the flux dilution, the geometry of 
the extended source is generally not apparent in the maps and 
the observer may not always be provided with a precise source 
geometry. Nevertheless, the signal to noise ratio does not de- 
crease rapidly with the integration region size once the source 
is fully contained in it. In order to check the global method on 
an astrophysical case, we analyzed the IBIS/ISGRI data on the 
Coma Cluster, a bright known cluster of galaxies with an X-ray 
size of order 30' diameter We applied this method, increasing 
the integration region until a maximum signal to noise ratio was 
reached. In that way, we were able to extract for the first time 
the global intensity of this extended source as well as interest- 
ing constrain ts on its morphology . The results are presented in 
a joint paper jRenaud et al. 20061) . In a forthcoming paper, this 
method will be also used on the IBIS/ISGRI observations on 
the Vela Junior SNR (Renaud et al. 2006, in preparation). 

Usually, the observer needs a lot of integration time for 
detecting an extended source (~ 5 x 10^ secondes for the 
Coma Cluster) and will have to construct mosaics of individual 
images. After obtaining the individual images with the standard 
analysis, the next step is to build mosaics as usual by weight- 
ing each image by their associated variance and transform them 
following the practical tips on the use of the method given in 
Appendix|B] Any observer interested in this method may con- 
tact the first author to obtain the map A. 
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Appendix A: Deconvolution Algorithm 

In this section, we will adopt the following notations: 

(k, I) detector 
(i, j) real sky 
(m, «) — > reconstructed sky 



(A.l) 
(A.2) 
(A.3) 



D is the detector, S and S' the real and reconstructed skies, 
M the coded mask, R the decoding matrix, and B the back- 
ground. By definition: 



Du = J^SijMi+t.j+, + B, 
ij 

S'„„ = Rk+mJ+nPkl 
kl 

with the following arrays: 



(A.4) 
(A.5) 



^k+m.l+n B^^nm X 



k+in,l+n 



k+mj+n 



ILkiG^ 



(A.6) 
(A.7) 



k+mJ+n 



'-'inn — - Bal„„, ^ G^_,„,_,_,„Ma:+,(,A,3) 

kl 



k+m,l+n 

kl 



where G^, G" the decoding arrays related to the coded mask M 
(see Goldwurm et al. 2003), Bal„,_„ is called the balance array 
in order to have a flat reconstructed flux image and E,„_„ corre- 
ponds to the imaging efficiency to correct for partial exposure 
outside the FCFOV. 



A.1. Reconstructed fluxes 

The standard reconstructed flux image is given by: 

^'mn - ^ Rk+mJ+n ^ S ijMi+kJ+l + ^ Rk+m,l+nBk,l (A. 9) 



We propose to construct images of flux per steradian (per 
sky pixel) I'mn = S'„„JA,„„ where A„,„ is the SPSF area. The 
global flux of any source (point-like or extended) I'^ is then 
given by the sum of the intensities within a region Q encom- 
passing the source extent. Therefore: 



4 - ^ -~ ^ S ij ^ Rk+mJ+nMi+kJ+l+y^^ -~ ^ Rk+m,l+nBkJ 

(A. 10) 

Assuming that the detector D was well background-substracted 
before the deconvolution, the second term of the right-hand 
side of the above equation is null. 

A„,„ represents the integrated normalized SPSF. By defini- 
tion, we have: 
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SPSFij{m,n) = ^Rk+mj+nMi+kj+i 

kl 

Amn = ^ SPSF,„„{m' ,n') 



(A.ll) 
(A. 12) 



where SPSF,j(m,n) is the SPSF value at the position (m,n) of 
the source at the position (ij) with a unit-flux. K'"" defines the 
local region encompassing the SPSF of the point source at the 
position (m,n). Substituting equation lA. 1 ll in lA. 101 



SPSFu{m,n) 



IJ m. 



ij m,neK'j 



SPSFij(m,n) 



(A.13) 
(A. 14) 



For each point source S,;, we have reduced the sum over 
Q (i.e. the source extent) to the sum over K'"-', the local region 
around the point source at (i,j). Thus, if we assume for each un- 
resolved point that A„„, = A,j is constant within K'-^, we finally 
obtain: 



(A. 15) 



A.2. Reconstructed variances 

From equation lA. 91 the reconstructed flux variance is given by: 



2 

k+m.l+n 



(A. 16) 

where cr^ . . is the variance related to the celestial source at the 
position (iJ) and cr^^., that related to the detector background. 
erf, ^jj^j is equivalent to the OSA flux variance o-'^^a- 

In most cases in the high-energy domain, the variance due 
to the celestial source(s) is very small compared to that due to 
the background count rate. Therefore, in the calcul of the vari- 
ance associated to the transformation I's - YjnS'„,„/A„,i,, the 
contribution of the first term in equation lA. lOl can be neglected. 
Given that the detector pixels are independent, we obtain: 



R 



k+m.l+n 



V n 



(A. 17) 



In the above equation, we have assumed for simplicity that the 

istant. As the deco 
1)+ and (2M - 1)-, re- 



detector variance cr^^^ is spatially constant. As the decoding 



arrays and G are given by (2M 
spectively, we obtain: 



kl ni,n2 



kl ^1,^2 



^ml,nl-^m2,n2Em2,n2 



(A.18) 



In the FCFOV, the balance array is almost -1 and only varies 
within ~ 5 % in the PCFOV. Therefore: 



kl ni.n. 



Rk+m\,l+n\ -Mk+mlJ+nl 
Am 1 ,nlAm2,n2F fii2,n2 
Rk+mlJ+nl-^ 



kl ili.n. 



A/iiljil Am2.n2F m2.n2 



(A. 19) 
(A.20) 



The second term of the equation is assumed to be 0, since 
the deconvolution of any flat image gives a reconstructed sky of 
values. We have expanded the first term, following the above 
definitions of A and R, and finally found a simple expression 
for the associated intensity variance such that: 



,2 2 Y' L 

I', ~ ^B,cl.s-^ Zj P a 

YjU Dk,i 



Q Finn A inn 



B,cl.s- 



P - 
J-'inn — 



Time^ 



On - Axis Aperture 



(A.21) 

(A.22) 
(A.23) 



In this expression, the final variance cr'? is in count s ^ for ho- 

mogeneity with I's (in count s"'). The On- Axis Aperture is the 
number of detector pixels which are illuminated through the 
mask by an on-axis source (a; half of the ISGRI camera, i.e. 
8192 pixels) and E^n is the imaging efliciency of the telescope. 

Appendix B: Practical tips on the use of the 
method 



Equation IA.21I presents a simple expression of the intensity 
variance associated with the method we developed for recon- 
structing the flux of an extended source. In order to provide the 
observer analyzing the IBIS data with an easy-to-use expres- 
sion, in relation to the standard OSA sky maps, we found that: 



B,ef..s-1 



and then, from equation lA.211 



OSA 



A mil 



(B.l) 



(B.2) 



If A does not suffer from large variations over Q, one may ap- 
proximate the last equation: 



' <A> 



■ N. 



spsf 



< a-. 



OSA 



(B.3) 



is an estimate of the OSA flux variance averaged 



where cr^, 

over Q.. The sum in equation lB.2l is reduced to the number of 
sky pixels within Q. divided by < A >, i.e. the number of IBIS 
SPSFs. The transformation 5 ^ /' allows the construction 
of sky images where anybody may integrate intensities over 
any region Q in order to properly measure the global flux of 
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an extended source. The associated variance is then given by 
equation IB. 21 The following equations allow the user to con- 
struct mosaics of a set of individual images: 



" mosaic n , / /9 



t^^Ka 



'^mosaic - ^ 1 , ,2 (^-^^ 

2j ^I"osa 



(B.5) 



As for an individual image, the global flux and the associated 
error of an extended source will be given by: 



= Y^IjI}£^ (B.7) 

Q ^mosaic 
,-'2 



n 



mosaic 



